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ABSTRACT 


An  experimental  study  of  the  radiation  properties  of  a 
parameter  fog  horn  was  conducted  using  a  4  unit  x  5  unit  horn 
dirver  array  as  the  primary  source.   The  horn  was  operated 
as  two  separate  10-unit  acoustic  arrays  each  transmitting 
one  primary  frequency.   The  frequency  combinations  studied 
were  4240-5520Hz  and  4960-5520HZ.   The  acoustic  properties 
of  the  parameter  array  were  observed  for  both  frequency 
combinations.   Source  levels  were  varied  between  122dB  re 
20iJPa  and  134dB  re  20yPa  at  1  meter  to  study  the  transition 
from  absorption  limited  performance  to  saturation  limited 
performance.   Results  consistent  with  the  predictions  of 
Moffett  and  Mellen  [JASA  61,  325-337  (3^977)]  for  the  source 
level  of  the  difference  frequency  beam  and  saturation 
broadening  of  the  difference  frequency  beam  were  observed. 
The  horn  was  also  operated  as  a  single  array  with  all  units 
radiating  both  primary  frequencies.   Equivalent  results  were 
obtained.   A  discussion  of  fog  signaling  applications  is 
included. 
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I.   INTRODUCTION 

A.   BACKGROUND 

The  U.S  Coast  Guard's  shore-based  sound-signal  system 
includes  over  600  signals,  many  of  which  are  located  near 
communities  and  residential  areas.   As  part  of  its  automation 
program  the  Coast  Guard  is  replacing  multif requency  air- 
driven  sound  signals  with  pure  tone  electric  air  oscillators. 
Installation  of  these  pure  tone  oscillators  has  increased 
noise  complaints  of  residents  near  the  signals.   One  signal, 
the  unbaffled  ELG300-02  (dual-emitter) ,    is  capable  of  pro- 
ducing levels  of  80dB(A)  as  far  as  one  half  mile  to  the  rear 
of  the  signal.   Environmental  Protection  Agency  guidelines 
state  that  60DB(A)  is  the  maximum  level  of  a  pure  tone  in 
urban  areas.  [1]   To  avoid  complains  levels  of  50dB(A)  or  less 
are  probably  desirable  in  suburban  and  rural  areas. 

Bolt,  Beneranek  And  Newman,  Inc.,  developed  and  improved 
baffle  for  the  ELG300-02  (dual  emitter) .   The  baffle  design 
goal  was  attenuation  of  30dB  to  40dB  to  the  rear  in  the  dark 
sector  of  the  signal.   A  baffle  design  featuring  combined 
acoustically  hard  and  soft  surfaces  was  developed  that  pro- 
vided 30dB  of  attenuation  in  the  dark  sector.  [2]   The  improved 
baffle  has  two  drawbacks:   Its  large  size  (10ft  x  18ft)  and 
its  cost  ($10,000  each). 


B.  OUTLINE  OF  INVESTIGATION 

The  purpose  of  this  investigation  was  to  test  a  model 
parametric  source.   A  parametric  source  could  ideally  provide 
an  acceptable  sound  level  in  the  bright  sector  of  the  signal 
with  little  or  no  noise  in  the  dark  sector  of  the  signal. 

Previous  work  on  theparametric  array  in  air  has  been 
limited.   Shealy  did  some  initial  work  on  parametric  beam 
widths  and  directivity  patterns.  [3]   Bennett  did  work  on 
propagation  of  the  difference  frequency,  however,  primary 
source  levels  were  not  high  (95dB  to  105dB  re  20yPa) ,  and 
the  difference  frequence  was  very  weak  (4  0dB  to  50dB  re 
20yPa) .  [4]   Brinkman  reported  a  difference  fequency  source 
level  of  108dB  re  20ijPa  from  primary  source  levels  of  130dB 
re  20]jPa.  [5]   Brinkman 's  difference  frequency  was  lOkHz 
which  was  too  high  for  fog  signaling  applications.   None  of 
the  previous  experiments  had  been  carried  out  to  the  farfield 
of  the  parametric  array.   Brinkman ' s  single  data  point  was  the 
only  one  available  to  check  the  Moffett  and  Mellen  design 
curves  in  air  [6]  and  there  is  some  question  of  how  close 
Brinkman 's  point  is  to  the  predicted  level.   It  was  concluded 
that  previous  work  on  the  parametric  array  in  air  was  not 
sufficient  to  predict  performance  in  fog  sinaling  applications 

C.  PERFORMANCE  REQUIREMENTS  FOR  FOG  SIGNALING 

A  set  of  minimum  performance  criteria  is  proposed  as  a 
measure  of  the  usefulness  of  the  parametric  array  in  air  for 
fog  signaling  applications: 


(1)  Difference  frequency  greater  than  lOOHz  but  less 
than  llOOHz  [7] 

(2)  Difference- frequency  level  greater  than  120dB  re 
20yPa  giving  a  "usual  range  of  detection"  of  the  signal  of 
one  mile  at  llOOHz  [8] 

(3)  Significant  reduction  of  noise  from  backward  radia- 
tion (The  baffle  design  goal  of  30dB  attenuation  in  the 
dark  sector  will  be  adopted.) 

(4)  Capability  to  steer  the  narrow  difference- frequency 
beam  over  large  angles. 


II.   PARAMETRIC  SOURCE  THEORY 

The  interaction  of  two  intense  sound  waves,  called 
primary  waves,  of  coincident  direction  of  propagation  and 
different  frequencies  results  in  the  formation  of  combina- 
tional waves.   The  combinational  waves  propagate  at  frequen- 
cies equal  to  the  sum  and  the  difference  of  the  primary 
wave  frequencies.   The  difference- frequency  wave,  also  called 
secondary  pressure  wave,  has  greater  engineering  application 
due  to  its  property  of  high  directionality  at  low  frequency. 

In  1963,  Westervelt  outlined  the  basic  solution  for  the 
secondary  pressure  wave  at  the  difference  frequency.  [9] 
The  secondary  pressure  wave  is  assumed  to  be  generated  by  a 
distributed  volume  of  sources  whose  strengths  depend  on  the 
local  level  of  the  primary  waves.  This  volume  source  is 
termed  the  "parametric  source."  The  spreading  geometry  and 
absorptive  behavior  of  the  primary  waves  in  the  interaction 
region  has  been  the  basis  for  different  theortical  develop- 
ments . 

A.   NEARFIELD  SOURCE        

Westervelt 's  original  theory  dealt  with  the  difference 
frequency  signal  generated  from  narrow,  plane,  collimated 
primary  waves.   It  is  the  model  for  a  parametric  source  con- 
strained by  absorption  of  the  primary  waves  in  the  nearfields 
of  the  primary  beams.   The  Rayleign  length,  R  ,  defines  the 
transition  range  from  the  nearfield  to  the  farfield  of  the 
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primary  beams.   The  Rayleigh  length  is  the  ratio  of  the 
source  area,  A  ,  to  the  mean  primary  wave  length,  A  . 


A  f 
o  o 


(2-1) 


where 


h^'i 


=  mean  primary  wave  frequency. 


Westervelt  found  a  solution  for  the  secondary  pressure 
wave 


B  CO,   P,  P^  A 
P,(r,0)   =  2 i 1 — o   exp(-a,  +  jk  )r  D(0)   (2-2) 

47Tr  p   c   a„ 
O      T 


where 


r,9  =  field  coordinates 

S   =  parameter  on  non-linearity  =  1.2  for  air 
a,   =  attenuation  coefficient  of  secondary  wave 
a^  ^  =  attenuation  coefficient  of  primary  waves 

=  o^  +  a2  -  a^ 

=  pressure  of  the  primary  waves 

=  density  of  medium 

=  speed  of  sound 

=   angular  frequency  of  secondary  wave 

27T 

=   T —  where  X-.    is  the  difference  frequency 
d 
wavelength. 


1.2 

1,2 

Po 
c 


and 
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2k   2     4  fl   -1/2 
DO)   =   (1   +   (— ^)^  sin^  f  )  (2-3) 


where  D(9)   is  the  directionality  factor  of  the  secondary 
wave. 

The  secondary  wave  is  generated  near  the  source  wherever 
the  product  of  the  primary  wave  pressures  P-iP^  is  large. 
Once  generated  the  secondary  wave  propagates  with  its  own 
attenuation  appropriate  to  that  frequency.   The  secondary 
sound  pressure  will  be  down  3dB  at  the  angle  6^  where 


.   ^d      ,  ^'t  .1/2  ,^  ., 

sm  -Y     =       (2]r~^  (2-4) 


if   ctm/k^   is  small,  the  beam  width  of  the  secondary  wave  is 


2  a 
d 


This  secondary  pressure  wave  is  a  narrow,  nearly  sidelobe- 
free  beam  of  low  frequency. 

B.   FARFIELD  SOURCE 

Berktay  and  Leahy  extended  Westervelt's  model  to  a 
farfield  source  where  the  bulk  of  the  interaction  is  assumed 
to  take  place  at  a  distance  greater  than  the  Rayleigh  length. 
[10]   In  this  region  the  primary  waves  are  assumed  to  spread 
spherically.   Berktay  and  Leahy  chose  to  write  the  on-axis 
secondary  pressure  in  terms  of  the  acoustic  power  W,  ^  °^  "^^^ 
primary  waves  where 
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"l,2   =  '^o'^,2>'/2PoC^  (2-6) 


The  secondary  pressure  becomes 


M  2w  1/2^  1/2   .(c^  +jk  )r 

P^Cr.O)   =   ^   ^  3   ^ e   ^   ^  (2-7) 

2ttc^  ra^ 


Substituting  the  appropriate  constants  for  air  and  writing 
a   in  terms  of  9,,  the  RMS  value  of  P,  at  Imeter, 

p    (1,0)       =       [_i41i80lLj    ^      ^      ,    ^  (2-8) 


From  this  expression  a  preliminary  design  equation  may  be 
deduced  as 


SL^  =   127.5  +  201ogf^  +  lOlogW^  +  lOlogW   -  401oge 

+  201og|v|  (dB  re  20yP  )  (2-9) 

3. 


where  f,  is  in  kHz,  W  ,W   are  in  watts,  0,  is  in  degrees  and 
201og|v|  is  a  correction  factor  involving  transducer  geometry. 
[10]   A  more  convenient  form  of  the  design  equation  results 
from  use  of  the  source  levels  of  the  primary  beams. 


SL.   =   lOlogW.  +  DI.  +  109 (dB  re  20yP^)         (2-10) 
1  11  ^  a 


In  air  this  yields 
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SL^   =   -90.5  +  201ogf,  +  (SL^-DI.)  +  (SL.+DI_) 
D  all        22 

-  401oge,  +  201og|v|(dB  re  20yP  )         (2-11) 

In  a  farfield  source  Westervelts  beam  width  is  only  good 
for  very  narrow  primary  beams.   If  20,  is  much  smaller  than 
the  product  of  the  primary  directivity  functions ,    then 
D(e)  becomes  the  product  of  the  primary  directivity  functions. 

C.   COMPREHENSIVE  MODEL 

There  are  two  additional  operating  regimes  for  parameter 
sources.   At  high  input  powers  the  parametric  array  length 
becomes  limited  by  nonlinear  absorption.   These  arrays  are 
called  saturation  limited,  because  their  effective  length  is 
determined  by  saturation  of  the  primary  signal  from  the 
generation  of  harmonics  and  shock  formation.   Moffett  and 
Mellen  have  developed  a  model  which  incorporates  all  four 
operating  regimes  of  the  parametric  array.  [6] 

The  on-axis  secondary  pressure  for  the  comprehensive 
model  is  written 

2     2 

ttBP^  R^f  , 

P^(r,0)   =   1— ^-^— (^1  +  ^2^  (2-12) 

PC  f^r 


where  the  integral  I^  is  a  measure  of  the  contribution  from 

the  collimated  portion  of  the  primary  wave  (Westervel  model) 

and  the  integral  I-  is  a  measure  of  the  contribution  from 

the  spherically  spreading  portion  (Beektay-Leahy  model) . 

The  source  area,  A  ,  of  Westervelt's  model  has  been  written 

o 
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in  terms  of  the  Rayleign  length,  R  ,  the  primaries  have  been 
assumed  to  have  the  same  pressure,  P^  =  P^  =  P.,  and 
frequency,  f^,  and  I,  and  I-  are  functions  of  ctm* 

The  complex  parameter  gain  is  defined  as  the  ratio  of  the 
on-axis  secondary  pressure  amplitude  to  the  amplitude  of  one 
primary 

g   E   rP^(r,0)/P^R^  (2-13) 


Substituting  P^  in  yields 


f   ^  In+In 

o      o 


where 


X   =   2Tr6PQRQfo/pc^  (2-15) 

In  decibel  form,  the  parametric  gain,  PG,  is  defined 

PG   =   20  log  |g|  (2-16) 

and  has  been  evaluated  and  plotted  on  a  series  of  design 
curves . 

The  parametric  gain,  PG,  is  a  function  of  three  dimen- 
sionless  parameters: 
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(1)  aR   .   The  product  of  the  mean  primary  wave 
absorption  (dB/m)  and  the  Rayleigh  length 

(2)  f.   /t^    .      The  ratio  of  the  difference  frequency,  f , , 
to  the  mean  primary  frequency,  f  .   It  is  called 
the  "downshift  ratio." 

(3)  X  .   In  lieu  of  x  it  is  often  more  convenient  to 
use  a  scaled  rms  primary  source  level,  L  *,  for 
each  primary  component  where 


P  R 
L  *   =   20  log  -^-^  +   20  log  f  (2-17) 

(kHz) (dB  re  20yPa-m-kHz) 


which  is  related  to  x  by 


20  log  X  =   L  *(dB  re  lyPa-m-kHz)  -  281  in  water 


(2-18) 


and 


20  log  X  =   L  *(dB  re  20yPa-m-kHz)  -  167  in  air 


(2-19) 


The  design  curves  of  Moffett  and  Mellen  can  thus  be  adopted 
for  air  by  replacing  281  (dB  re  lyPa-m-kHz)  with 
167  (dB  re  20yPa-m-kHz)  [15]. 

In  using  the  design  curves,  the  difference  frequency 
source  level  is  found  by  adding  the  parametric  gain  to  the 
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actual  input  source  level  utilizing  the  appropriate  downshift 
ratio,  scaled  input  source  level,  and  absorption  parameter. 
The  directivity  gain,  DG,  from  the  curves  is  added  to  the 
directivity  index  of  the  primary  waves ,  DI  ,  to  obtain  the 
directivity  index  of  the  secondary  wave,  DI^. 

DI,   =   DI   +  DG  (2-20) 

d       p 

As  an  example  of  how  the  curves  are  utilized,  consider  the 
following  system.   A  source  driven  at  primary  frequencies 

4.960kHz  and  5.52kHz.   We  have   f   =  ^-96  +  5.52  ^  5.24kHz, 

o         2 

f,  =  0.56kHz  and  downshift  ratio,  f  /f ,  =  9.3.   If  the  area 
d  o^  d 

2 
of  the  source  =  0.26m.  ,  then  the  Rayleigh  length,  R  =4.02. 

-2  -3 

For  frequency,  f  ,  we  have  a  =  3.1x10   dB/m  =  3.56x10   np/m 

[16]  and  otR  =  0.12dB.   If  the  source  level  of  each  of  the 

primaries  is  130dB  re  20yPa,  then 

L  *   =   scaled  primary  source  level 

P  R 
=   20  log  -2-^  +  20  log  f   (kHz) 
/2"  ° 

=   130  +  20  log  R  +  20  log  5.24 
=   156.4  (dB  re  20yPa-m-kHz) . 

Now,  use  of  Figure  2-1  yields  parametric  gain,  PG  =  -37.5dB 
and  the  predicted  secondary  level  is 
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p 

SL,   =   20  log  -——  +   PG 
a  P   JT 

ref 


=   130  -  37.5 


=   92.5  (dB  re  20yPa) 


Assuming  the  same  directivity  index  for  each  of  the 
primaries,  from  Figure  2-1,   DG  =  0   and   DI,  equals  DI 
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°  DIRECTIVITY 


117  lyr  14-7  157  167  177 

SCALED  PRIMARY  SOURCE  LEVEL 
{dBll20fi  PA-rn-  kHii) 

FIGURE  2-1.   MOFFETT  AND  MELLEN  DESIGN  CURVE 
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III.   EXPERIMENTAL  PARAMETRIC  SOURCE 

A  source  capable  of  high  output  over  a  broadband  of  fre- 
quencies was  considered  necessary  to  carry  out  a  useful 
parametric  experiment.   Available  at  the  Electroacoustics 
Research  Laboratory,  University  of  Texas,  Austin  and  owned 
by  the  U.S.  Coast  Guard  was  an  array,  manufactured  by  Applied 
Electro  Mechanics  (AEM) ,  Alexandria,  Va.   This  AEM  array 
had  been  extensively  tested  at  the  University  of  Texas, 
Austin  and  the  basic  linear  parameters  were  reported  in 
reference  [11]  and  reference  [12] .   They  are  reproduced  here. 

A.   ARRAY  CHARACTERISTICS 

The  array  is  rectangular  and  composed  of  20  individual 
exponential  horns.   Each  horn  flares  exponentially  from  a 
round  throat  2.5cm  in  diameter  to  a  square  mouth  11.4cm 
on  a  side.   The  horn  length  is  35cm  and  each  horn  is  driven 
by  a  separate  35watt  horn  driver.   The  horns  are  bolted 
together  to  form  a  four  element  by  five  element  array  with 
dimensions  45.6cm  by  57.5cm. 

The  relative  frequency  response  of  the  array  is  shown  in 
Figure  3-1.   The  array  responds  from  250Hz  to  5800Hz  with 
peaks  at  1500Hz,  3500Hz,  and  5300Hz.   The  beam  patterns  of 
the  array  at  lOOOHz  and  2000Hz  are  shown  in  Figure  3-2.  [11] 
At  2000Hz,  the  array  has  become  very  directional  with  very 
little  backward  radiation.   Unbaffled  it  very  nearly  meets 
the  performance  requirements  of  30dB  attenuation  in  the  dark 
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sector.   A  similar  lack  of  backward  radiation  can  be  expected 

at  frequencies  higher  than  2000Hz.   Figures  3-3,  3-4,  and 

3-5  are  the  main  lobe  beam  patterns  across  the  4  element  side 

at  the  frequencies  of  4240Hz,  4960Hz,  and  5520Hz.   These 

frequencies  will  be  used  to  generate  the  parametric  beam. 

All  three  beam  patterns  are  symmetric  with  -3dB  beam  widths 

within  experimental  error  of  theoretical  predictions. 

Reference  12  reported  that  the  array  was  capable  of 

producing  sound  pressure  levels  at  one  meter  of  14  6dB  re 

20ijP^  at  1500Hz,  144dB  re  20yP,  at  3500Hz,  and  142. 5dB  re 
a  a 

20yP   at  5333Hz.   These  sound  pressure  levels  were  produced 
from  an  input  of  approximately  600watts.   The  high  source 
levels,  broadband  response,  and  reduced  backward  radiation 
above  lOOOHz  made  the  array  an  attractive  model  parametric 
source. 

B.   PARAMETRIC  OPERATION  OF  THE  ARRAY 

The  array  was  employed  in  two  different  ways  to  get  an 
acceptable  region  of  interaction  for  parametric  generation 
of  sound. 

1.   Mixed  Method 

Both  frequencies  were  radiated  from  all  twenty 
drivers  of  the  array.   By  setting  each  oscillator  at  a  primairy 
frequency  and  using  a  resistance  bridge  to  mix  the  signals, 
a  combined  signal  was  produced.   The  signal  was  high  passed 
through  a  filter  to  keep  only  the  primaries,  then  amplified, 
and  sent  to  the  array  (Fig.  3-6) .   Care  had  to  be  taken  to 
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FIGURE  3-4.   4960H2  PRIMARY  BEAM  PATTERN 


FIGURE  3-5.   5520Hz  PRIMARY  BEAM  PATTERN 
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match  the  output  impedance  of  the  amplifier  to  the  input 
impedance  of  the  array  for  two  reasons:   (1)  impedance 
matching  maximizes  the  amount  of  power  transferred,   (2)  it 
was  noted  that  a  failure  to  impedance  match  resulted  in  the 
horns  directly  radiating  the  difference  frequency  through 
intermodulation  distortion. 
2.   Split  Method 

In  this  method  two  oscillators  (each  set  at  a  primary 
frequency),  two  amplifiers,  and  the  acoustic  sources  were 
utilized.   The  array  functioned  as  two  separate  sources  by 
wiring  ten  drivers  to  one  amplifier  and  ten  drivers  to  the 
other  amplifier.   Each  source  formed  a  checkerboard  pattern 
on  the  array  face.   Figure  3-6.  The  primary  waves  were 
electronically  isolated  and  only  mixed  as  acoustic  waves  in 
the  air. 


C.   TRANSMITTING  EQUIPMENT 

A  list  of  the  transmitting  equipment  utilized  is  given 
below: 


(1)   Oscillators 


(2) 
(3) 


Filter 

Power 
Amplifier 


(4)   Power 

Amplifier 


2  ea.  GR  type  1163-A  decade 
frequency  synthesizer 

Krohn-Hite  Model  3202R  filter 

Kilowatt  amplifier  model  LDV2-3 
watt  with  selectable  matching 
load  resistances 

Kilowatt  amplifier  model  LDC3-2 
520  watt  with  selectable  matching 
load  resistances. 
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The  kilowatt  amplifiers  had  several  capabilities  which 
were  especially  useful  for  these  experiments.   For  an  input 
of  IVolt  RMS  they  will  produce  full  rated  power  output, 
eliminating  any  preamps  which  could  introduce  distortion. 
A  1:100  voltage  divider  across  the  output  terminals  of  the 
amplifier  is  brought  to  a  BNC  voltage  monitor  connector  on 
the  front  face  of  the  amplifier.   The  voltage  from  the 
monitor  is  low  enough  to  be  spectrum  analyzed  to  test  output 
signal  purity.   The  array  was  wired  for  testing  as  two 
circuits  of  ten  drivers  each.   All  ten  drivers  in  each 
circuit  were  in  series  to  raise  the  impedance  of  the  network 
to  match  the  output  impedance  of  the  amplifier.   By  wiring 
the  drivers  in  series  the  failure  of  any  driver  would  break 
the  circuit  giving  an  easily  detectable  warning  of  degraded 
array  performance. 
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IV.   RECEPTION  AND  RECEIVING  EQUIPMENT 

A.   RECEPTION  OF  PARAMETRIC  SIGNALS 

In  receiving  parametric  signals  it  is  necessary  to  insure 
that  extraneous  difference  frequency  signals  do  not  interfere 
with  the  actual  parametrically  generated  difference  frequency 
signal.   Bennett  identified  three  possible  sources  of 
extraneous  signals  and  their  properties  [4].  They  are: 

(1)  Intermodulation  distortion  in  the  receiving 
equipment.   This  extraneous  signal  varies  as  the  product  of 
the  local  primary  amplitudes  at  the  microphone  and  exhibits 

a  propagation  curve  whose  slope  is  12dB  per  doubling  of  range, 

(2)  Pseudo  sound  interaction.   This  is  a  signal  produced 
on  the  face  of  the  receiver  by  radiation  pressure.   It  also 
exhibits  a  slope  of  12dB  per  doubling  of  range. 

(3)  Transmission  distortion.   This  signal  can  be 
identified  by  its  broad  beam  directionality.   By  utilizing 
the  split  transmitting  method  this  source  of  extraneous 
signal  is  eliminated. 

The  behavior  of  the  spurious  signals  is  contrasted  with 
the  parametrically  generated  difference  frequency  whose 
propagation  curve  first  rises  with  distance,  reaches  a  broad 
maximum,  and  then  gradually  decreases  approaching  a  slope 
of  6dB  per  doubling  of  range  as  the  farfield  of  the  array 
is  reached  [4] . 
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B.   RECEIVING  EQUIPMENT 

Three  (different  receiving  systems  were  tested.   The  goal 
was  to  finci  the  system  that  minimized  the  extraneous  difference 
frequency  signals.   Characteristics  of  the  three  systems 
tested  are  listed  below. 


Receiver 

1/2"  ALTEC 

BR-150 

Microphone 

1/2"  ALTEC 
BR-200  Mi 

Sensing 
type 

Condenser 

Condenser 

Sensitivity 
re  1  volt/ 
ybar 

-58dB 

-85dB 

Linear 
Limit 

158dB 

185dB 

Average  Noise    58dB 

Threshold 

SPL  re  20yPa 

85dB 

Frequency 
Response 

lOHz  - 
6500Hz 

30Hz- 
6000Hz 

Power 
Supply 

ALTEC 
526B 

ALTEC 
526B 

Amplifier 

none 

none 

LC-10  Hydrophone 


Lead  Zirconate 
Titanate 

-108dB 


iHz  to  12  0kHz 


none 


HP465A  with 

40dB 

amplification 


Each  receiving  system  was  configured  identically  to  its 
proposed  method  of  operation  during  actual  parametric 
testing.   For  the  ALTEC  BR-150  and  BR-200  the  system  included 
the  microphone,  power  supply  and  150'  of  microphone  cable. 
The  LC-10  system  included  the  hydrophone,  amplifier  and 
150'  of  cable. 
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To  measure  the  intermodulation  of  the  microphone  systems 
the  array  was  configured  in  the  split  transmitting  method. 
This  eliminated  the  chance  of  any  transmitted  spurious 
difference  frequency  signal.   Each  microphone  was  placed 
less  than  0.5m  from  the  array.   The  array  was  turned  on, 
and  the  levels  of  the  primaries  received  by  the  microphones 
were  varied.   As  the  microphones  were  very  close  to  the 
source  it  was  assumed  that  any  difference  frequency  being 
recorded  was  the  result  of  intermodulation  in  the  receiving 
system  and  not  parametrically  generated.   This  was  a  safe 
assumption  as  the  overlap  distance  of  the  two  adjacent 
units  in  the  split  system  was  .41m.   The  spurious  difference 
frequency  generation  in  the  receiving  equipment  and  pseudo- 
sound  could  be  measured  in  total.   The  results  are  plotted 
as  equivalent  sound  pressure  level  of  the  spurious  signal 
vs.  sound  pressure  level  of  the  primary  signals  at  the 
microphone.   Fig.  4-1.   Also  noted  is  the  noise  threshold 
of  each  system. 

The  LC-10  hydrophone  was  chosen  as  the  receiving  system. 
For  a  given  primary  input  it  had  30dB  less  intermodulation 
than  the  BR-150  and  5dB  less  than  the  BR-200.   The  LC-10 
observed  noise  threshold  was  13dB  less  than  the  BR-150  and 
40dB  less  than  the  BR-200.   A  primary  signal  to  difference 
frequency  intermodulation  ratio  of  60dB  established  the 
lower  limit  of  the  level  of  parametric  signals  which  could 
be  measured  by  the  LC-10. 
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The  remainder  of  the  receiving  equipment  included: 

(1)  Spectrum  analyzer    Schlumberger  Model  1510-03 

(2)  Wave  analyzer        HP  Model  30 2A 

(3)  Filter  Krohn  Kite  Model  3202R 

(4)  Sound  level  GR  Type  1562-A 
calibrator 

The  Schlumberger  Model  1500-03  was  utilized  as  the  level 
recording  device.   The  model  1510-03  is  a  fully-digitalized 
instrument  that  provides  real  time  spectral  analysis  of 
analog  signals.   It  has  five  full  scale  input  ranges  from 
O.lV  full  scale  to  lOV  full  scale.   The  input  signal  to  the 
model  1510-03  is  divided  into  256  spectral  bands  (lines  of 
resolution)  across  one  of  10  selectable  frequency  bands  of 
interest.   The  frequency  spacing  of  the  lines  of  resolution 
vary  from  O.lHz  to  lOOHz.   For  one  testing  on  the  primary 
frequencies  (5000Hz)  a  40Hz  spectral  bandwidth  with  a  12.8s 
averaging  time  was  used.   For  the  difference  components,  at 
1280Hz,  a  lOHz  spectral  bandwidth  with  6.4s  averaging  time 
and,  at  560Hz,  a  4Hz  spectral  bandwidth  with  16s  averaging 
time  were  utilized.   To  measure  the  correct  levels  care  was 
taken  to  choose  frequencies  that  placed  the  primary  and 
secondary  frequencies  in  the  center  of  the  rectangular 
sampling  window  of  the  spectriom  analyzer.   The  difference 
frequencies  were  also  chosen  to  coincide  with  spectral  bands 
that  had  low  noise  levels. 
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C.   DATA  ACQUISITION 

The  initial  testing  of  the  array  at  the  Naval  Postgraduate 
School  was  done  in  the  anechoic  chamber.   The  maximum  distance 
that  a  signal  could  be  propagated  in  the  chamber  was  7m. 
The  far field  of  the  array  is  generally  considered  to  begin 
at  l/2a  where  a   =  average  primary  absorption  in  np/m. 
For  a  5000Hz  signal  the  distance  to  the  farfield  is 
approximately  125m.   The  7m  propagation  distance  was 
considered  to  be  too  short  to  get  useful  data  so  the 
experiment  was  moved  outdoors.   The  outdoor  experiments 
i  were  done  using  the  Mars  Radio  Tower  on  the  grounds  of  the 
Naval  Postgraduate  School.   The  array  was  placed  facing 


upward  2.5m  from  the  tower  on  a  level  stand.   The  stand  was 
leveled  prior  to  each  data  run  with  a  carpenter's  level 
placed  on  the  face  of  the  array.   The  level  was  measured  to 
be  accurate  to  within  1/4*'.   The  receiving  microphone  was 
attached  to  a  boom  and  extended  2 . 5m  out  from  the  tower .   A 
system  of  two  surveying  theodelites  was  used  to  position 
the  microphone  directly  over  the  center  of  the  array.   The 
positional  accuracy  of  the  microphone  was  within  2cm  of 
the  acoustic  axis  of  the  array. 

The  propagation  distance  was  varied  by  moving  the  boom 
to  different  heights  on  the  tower.   There  were  only  four 
possible  propagation  distances  the  maximum  of  which  was  30m. 
Scattering  from  the  tower  was  ignored  as  the  dimensions  of 
the  tower  stock  were  much  less  than  the  wavelengths  of  the 
primaries . 
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To  take  beam  patterns  outdoors  the  array  was  placed 
on  a  tiltable  platform.   It  was  leveled  using  a  carpenter's 
level.   To  measure  angles  a  12-inch  protractor  was  taped 
to  the  side  of  the  array  with  the  protractor's  base  level 
with  and  parallel  to  the  face  of  the  array.   A  line  and 

eI  plumbob,  attached  at  the  apex  of  the  protractor  and 

calibrated  to  read  0°  when  the  face  of  the  array  was  level, 
were  used  to  measure  angles  of  tilt.   The  maximum  angle 

\   that  could  be  measured  with  this  system  was  about  5°   . 

For  more  complete  beam  patterns  the  array  was  placed  in  the 

I    anechoic  chamber  at  the  Naval  Postgraduate  School.   A  small 
motor  turns  the  source  while  simultaneously  sending  a  signal 
to  the  X-scale  on  the  plotter.   The  level  of  the  signal  from 
a  microphone  mounted  in  the  chamber  is  sent  to  the  Y-scale 
of  the  plotter. 

D.   EQUIPMENT  CHECKOUT  AND  CALIBRATION 

Prior  to  the  collection  of  data  the  following  checks 
were  made: 

(1)  The  LC-10  hydrophone  was  comparison  calibrated 
with  a  BR-150  microphone.   The  BR-150  was_fir3t  calibrated 
with  a  GR-1652-A  piston  calibrator.   The  LC-10  and  the  BR-150 
were  then  placed  at  the  same  position  in  a  sound  field  and 
the  responses  compared. 

(2)  The  Schlumberger  spectrum  analyzer  was  calibrated 
using  its  internal  signal  generator. 
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(3)  The  array  was  precisely  positioned  using  the 
theodelites  and  the  carpenter's  level. 

(4)  The  array  was  checked  for  driver  failure  by  simply 
sending  a  signal  through  each  circuit  and  listening  for  a 
response. 
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V.   EXPERIMENTAL  INVESTIGATION 

Initially,  the  purpose  of  this  investigation  was  to 
test  the  array  in  a  fog  signaling  configuration  where  the 
primaries  would  be  set  around  2000Hz  and  the  difference 
frequency  would  be  around  400Hz.   This  configuration  had 
to  be  abandoned  as  the  testing  generated  noise  complaints 
from  the  occupants  of  residential  areas  near  the  testing 
tower.   In  order  to  continue  testing  without  causing  more 
noise  complaints/  the  primary  frequencies  were  raised  to 
the  vicinity  of  5000Hz.   The  noise  complaints  stopped  after 
this  step  was  taken. 

It  was  concluded  that  this  research  would  concentrate 
on  the  verification  of  the  ability  of  the  Moffett  and  Mellen 
design  curves  to  predict  the  source  level  of  the  difference 
frequency  in  air.   Only  the  single  data  point  from  Brinkman ' s 
experiment  was  previously  available  to  compare  with  the 
Moffett  and  Mellen  predictions. 

A.   PROPAGATION 

The  Moffett  and  Mellen  design  curves  are  farfield 
predictions  for  the  parametric  array.   Based  on  a  5000Hz 
primary  frequency,  a  farfield  distance  of  130m  was  expected 
for  the  experiment.   As  the  maximum  range  of  observation 
was  30m,  all  measurements  were  limited  to  the  nearfield  of 
the  parametric  array.   Theoretically,  the  difference  frequency 
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propagation  curve  in  the  nearfield  should  first  rise  with 
distance,  reach  a  broad  maximum  and  then  decrease  approaching 
a  slope  of  6dB  per  doubling  of  range  as  the  farfield  of  the 
array  is  reached.   Mellen  [13]  and  Rolleigh  [14]  have 
developed  correction  factors  for  the  nearfield.   Plotted 
in  Figures  5-1  and  5-2  are  the  extrapolated  farfield 
spreading  curves  as  predicted  by  Moffett  and  Mellen.   Below 
this  curve  is  the  Rolleigh  nearfield  propagation  curve. 
Rolleigh 's  result  is  limited  to  the  spherical  spreading 
region  of  the  primaries,  which  is  greater  than  4m  in  both 
cases.   Mellen 's  prediction,  although  not  plotted,  gives  a 
similar  curve. 

The  data  plotted  in  Figures  5-1  and  5-2  are  observed 
from  the  split  transmitting  method.   At  short  range  the  data 
appear  below  the  predicted  level  but  at  longer  ranges  they 
approach  the  Rolleigh  prediction.   This  behavior  is  probably 
the  result  of  the  split  transmitting  method.   The  beam  width 
of  each  individual  horn  is  33°  at  5000Hz.   Not  until  the 
sound  has  propagated  approximately  2.0m  do  the  beams  of  all 
units  overlap.   This  would  indicate  that  the  first  2.0m  of 
the  interaction  region  is  somewhat  weakened.   This  loss  of 
interaction  region  is  not  important  at  longer  ranges  because 
this  is  a  farfield  source  where  the  bulk  of  the  interaction 
takes  place  in  the  farfield  of  the  primary  waves.   The  loss 
of  the  first  of  the  interaction  region  in  the  nearfield 
could  be  responsible  for  the  low  readings  at  short  ranges. 
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B.   SOURCE  LEVEL  DIFFERENCE  VS.  SOURCE  LEVEL  PRIMARY 

In  reviewing  the  Moffett  and  Mellen  design  curves  it 
can  be  noted  that  at  low  primary  source  levels  a  3dB  increase 
in  the  primaries  source  level  produces  a  3dB  increase  in  the 
parametric  gain  or  a  6dB  increase  in  the  secondary  source 
level.   This  behavior  is  predicted  by  Westervelt,  and  Berktay 
and  Leahy,  for  an  absorption  limited  parametric  array  (where 
the  level  of  the  secondary  varies  as  the  product  of  the 
primary  levels) . 

At  high  levels  of  3dB  increase  in  the  primaries  produces 
less  than  a  6dB  increase  in  the  secondary.   This  is  the 
saturation  region  of  the  array.   In  this  region  energy 
instead  of  being  coupled  into  the  difference  frequency,  as 
is  done  at  low  levels,  is  being  lost  in  the  generation  of 
harmonics  of  the  primaries. 

The  array  had  enough  flexibility  in  source  level  to 
compare  to  the  saturation  effects  with  the  predictions  of 
Moffett  and  Mellen.   To  observe  this  behavior  the  microphone 
was  placed  on  the  propagation  axis  at  27.5m  and  the  output 
of  the  primaries  varied  from  low  to  high  levels.   The 
Moffett  and  Mellen  source  level  predictions  are  for  the 
farfield  region  while  the  observation  point  at  27.5m  was  in 
the  nearfield.   +4dB  was  added  to  the  data  points  at  12  80Hz 
and  +5dB  was  added  to  the  points  at  560Hz.   These  are  the 
corrections  predicted  by  Rolleigh  for  the  nearfield  at  27.5m. 
With  this  correction  and  extrapolating  both  the  difference 
frequency  and  the  primaries  back  to  one  meter,  it  is  possible 
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to  make  a  plot  of  source  level  difference  frequency  versus 
source  level  of  one  primary.   See  Figures  5-3  and  5-4. 
Illustrated  on  the  graph  are  the  Moffett  and  Mellen 
predictions  and  a  prediction  for  the  level  in  a  strictly 
absorption  limited  case. 

To  verify  that  saturation  is  taking  place  pictures 
were  taken  of  the  spectrum  at  6.5m  for  low  primary  source 
level  of  122dB  re  20yPa  and  high  primary  source  level  of 
130dB  re  20iJPa.   The  top  photograph  in  Figure  5-5  is  the 
spectrum  recorded  at  6.5m  from  the  low  primary  source  level, 
the  bottom  photograph  is  the  spectrum  from  the  high  level 
input.   Visable  to  the  right  of  the  two  primaries  in  each 
photograph  are  the  second  harmonics  and  the  sum  frequency. 
At  the  higher  primary  levels  the  second  harmonics  are  7dB 
closer  to  the  first  harmonic  levels  than  at  the  low  level, 
which  represents  a  loss  of  energy  from  the  interaction  region, 

By  shortening  the  interaction  region,  saturation  also 
broadens  the  difference  frequency  beam  patterns.   Figure  5-5 
illustrates  the  amount  of  broadening  that  occurred  in  the 
difference  frequency  beam  pattern  between  low  level,  120dB 
re  20yPa,  and  high  level,  130dB  re  20yPa.   At  low  levels, 
the  difference  frequency  beam  pattern  should  be  the  product 
of  the  primary  beam  patterns  [10].   At  high  levels,  the 
difference  frequency  beam  pattern  is  broader  and  can  be 
estimated  from  the  Moffett  and  Mellen  curves  [6] . 
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FIGURE  5-5. 


LOW  SOURCE  LEVEL  AND  HIGH  SOURCE  LEVEL 
SPECTRUMS  AT  6.5m 
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The  theoretical  pattern  for  the  low  level  was  calculated 
by  multiplying  the  theoretical  directivities  of  a  4240Hz  beam 
j  and  a  5520Hz  beam  calculated  from  the  dimensions  of  the 
5  unit  side  of  the  array.   The  theoretical  pattern  for  the 
high  level  was  found  assuming  OdB  directivity  gain  from  the 
design  curves  and  adding  it  to  an  average  directivity  of 
the  two  primaries.   Both  precictions  are  consistent  with  the 
data.   Figure  5-7  is  a  wider  view  of  the  1280Hz  difference 
frequency  beam  taken  at  6.5m  which  illustrates  the  narrow, 
sidelope-free  beam  shape  of  the  difference  frequency. 

C.   SPLIT  AND  MIXED  TRANSMITTING 

One  concern  of  this  experiment  was  determining  the 
equivalency  of  the  split  and  mixed  transmitting  systems. 
Two  experiments  had  failed  to  get  an  acceptable  region  of 
interaction  from  a  horn  driver  system  in  air  configured  as 
a  split  transmitting  system  [3] ,  [4] .   As  noted  in  the 
previous  section,  it  appeared  that  the  first  2m  of  the 
interaction  region  was  being  lost  with  the  split  system. 
This  loss  was  found  not  to  be  significant  when  the  source 
level  primary  versus  source  level  secondary  data  runs  were 
made  again  with  the  array  in  a  mixed  transmitting  system 
(Figures  5-3,  5-4).   In  comparing  the  two  systems,  no 
trend  is  evident  in  the  data  to  suggest  an  effect  from  split 
transmitting  when  the  data  are  taken  at  27.5m. 
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FIGURE  5.6.   12  80Hz  DIFFERENCE  FREQUENCY  BEAM  at  2  7.5m 
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FIGURE  5-7.   12  8  0Hz  DIFFERENCE  FREQUENCY  BEAM  AT  6.5m 
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VI.   RESULTS  AND  CONCLUSIONS 

A.   EXPERIMENTAL  PROPERTIES  AND  THEORETICAL  PREDICTIONS 

This  experimental  work  on  the  parametric  array  in  air 
was  able  to  demonstrate  several  properties  of  the  array  that 
had  not  been  verified  before: 

(1)  An  acceptable  region  of  interaction  from  two 
electrically  isolated  sources  of  horn  drivers  generated  a 
difference  frequency  which  behaved  almost  identically  to  an 
equivalent  mixed  electronic  system  (Figs.  5-3,  5-4). 

(2)  The  slope  change  in  the  source  level  curves  of 
Moffett  and  Mellen  was  observed  when  the  source  bagan  to 
saturate.   By  applying  the  nearfield  correction  of  Rolleigh, 
excellent  agreement  was  shown  between  the  Moffett  and  Mellen 
predicted  levels  and  the  experimental  data  (Figs.  5-3,  5-4). 

(3)  The  behavior  of  the  parametric  array  in  the  nearfield 
was  found  to  follow  the  nearfield  predictions  of  Rolleigh. 

At  the  shortest  ranges  the  data  are  below  the  predicted 
level.   This  is  probably  the  result  of  using  the  mixed 
transmitting  system  where  the  interaction  region  at  distances 
of  less  than  2m  are  lost.   This  loss  does  not  affect  the 
overall  beahvior  of  the  array  as  the  bulk  of  the  interaction 
takes  place  at  a  distance  greater  than  4m,  but  could  effect 
the  propagation  curve  in  that  region  of  the  nearfield  close 
to  the  source  (Figs.  5-1,  5-2) . 


50 


(4)   Difference  frequency  beam  broadening  was  observed 
as  the  array  length  became  limited  by  saturation  (Fig,  5-5) . 

B.   FOG  SIGNALING  CAPABILITY 

Although  it  was  not  possible  to  test  the  array  in  an 
optimal  fog  signaling  configuration  the  Moffett,  Mellen 
design  curves  make  it  possible  to  propose  a  field  configura- 
tion for  a  parametric  fog  horn. 

2 
The  proposed  horn  would  be  a  50-unit  array,  area  0.5m  , 

operating  with  primaries  around  2000Hz  and  secondary 

around  500Hz,  input  power  1200  watts,  and  source  level 

primaries  of  145dB.   Entering  the  design  curves: 

Primary  =  2kHz    Secondary  =  50  0Hz    Downshift  =  4 
«2000  =  0-009dB/m      R^  =  -""340°°   =  2.9m 

aR^  =  0.02dB 

Scaled  Source  Level  =  145  +  20  log  2.9  +  20  log  2 

=  160dB  re  20yPa-m-kHz 

From  Fig.  Z-1   ,   parametric  gain  =  -20dB 

Difference  Frequency  Source  Level  =  125dB  re  20ijPa. 

This  proposed  parametric  fog  horn  will  be  compared  with 
a  linear  horn  operating  at  500Hz.   In  general,  the  amount  of 
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unwanted  sound  in  the  dark  sector  of  the  fog  signal  will 
depend  on  several  factors.   These  include  source  directionality, 
baffle  attenuation,  atmospheric  absorption,  geometric  spreading, 
physio  acoustic  effects  of  the  ear  receiver  and  psycho 
acoustic  perception. 

The  difference  frequency  beam  is  narrow,  nearly  side- lobe 
free.   Any  noise  from  the  parametric  fog  horn  will  be  from 
the  primaries  at  2000Hz,  while  noise  from  the  linear  source 
will  be  at  500Hz.   To  get  the  same  source  level  at  500Hz 
the  parametric  primaries  will  have  to  be  20dB  higher  than 
the  linear  signal.   At  the  face  of  the  array  the  parametric 
fog  horn  will  be  20dB  higher  in  level. 

The  directionality  of  the  array  reduces  the  parametric 
signal  behind  the  array.   The  array  is  nearly  omni- 
directional at  500Hz  but  at  2000Hz  it  is  30dB  down  in  the 
backward  direction. 

Both  sources  could  be  baffled  to  further  reduce  backward 
radiation.   The  parametric  signal  due  to  its  shorter  wavelength 
would  require  a  smaller  baffle  than  the  linear  signal  for 
similar  attenuation. 

Both  signals  would  spread  geometrically  at  the  same  rate, 
but  the  parametric  signal  would  be  attenuated  faster  due  to 
its  greater  atmospheric  absorption.   The  atmospher  absorbs 
500Hz  at  26dB/mile  and  2000Hz  at  the  rate  of  15.9dB/mile. 

The  physioacoustic  and  psychoacoustic  effects  are  the 
most  difficult  to  predict.   It  is  known  that  the  human  ear 
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is  more  sensitive  at  2000Hz  than  at  500Hz  giving  a  perceived 
sound  level  gain  to  the  parametric  signal.   Conversely  the 
parametric  signal,  which  is  made  up  of  two  frequency 
components,  with  proper  frequency  selection  could  present 
a  much  more  pleasing  pitch  raising  the  level  of  tolerance 
of  the  signal  above  a  simple  pure  tone  of  equal  intensity. 

The  parametrically  generated  difference  frequency  beam 
is  quite  narrow  compared  to  the  angular  width  a  fog  signal 
normally  illuminates.   To  compensate  for  the  narrow  beam 
the  source  would  have  to  be  scanned.   This  can  be  accomplished 
by  mechanically  moving  the  entire  array  assembly  through  the 
angles  desired.   Another  method  would  be  to  electronically 
scan  the  signal. 

In  the  experimental  work  all  horn  drivers  were  driven  in 
phase  which  creates  a  beam  pattern  where  the  major  lobe  is 
perpendicular  to  the  length  of  the  array.   By  introducing 
a  time  delay,  i.e.,  phase  difference  between  the  acoustic 
signals  generated  at  neighboring  drivers,  the  major  lobe 
can  be  pointed  in  an  arbitrary  direction.   With  a  micro- 
processor control  unit  the  time  delays  could  be  continuously 
varied  to  steer  the  beam  over  any  angle  desired. 

One  capability  from  narrow  beams  and  steering  not 
available  from  the  present  omnidirectional  fog  horns  is  the 
ability  to  vary  the  signal  content  to  different  sectors. 
For  example  a  signal  could  be  programmed  to  transmit  a 
short,  then  long  blast  to  sectors  to  the  left  of  a  naviga- 
tion channel,  a  long  blast  right  down  the  channel  and  a  long 
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then  short  blast  to  sectors  to  the  right  of  a  channel. 
Similarly,  a  shoal  could  have  a  different  signal  transmitted 
over  it  than  the  rest  of  the  sectors,  as  is  done  with 
sectored  lighthouses. 

C.  DISADVANTAGES  OF  A  PARAMETRIC  FOG  HORN 

The  primary  disadvantage  of  a  parametric  fog  horn 
would  be  cost  of  acquisition  and  operation.   The  proposed 
fifty  horn  unit  with  associated  amplifiers  and  oscillators 
could  cost  $25,000  each.   Parametric  generation  of  sound  is 
also  very  inefficient.   To  generate  a  125dB  re  20yPa  signal 
linearly  would  take  less  than  ten  watts  input  whereas  the 
parametric  input  would  be  1200  watts.   Additionally,  due  to 
the  need  to  scan  the  narrow  beam  over  large  angles ,  the  signal 
would  have  to  be  transmitting  roughly  ten  times  as  long 
during  each  broadcast  cycle. 

D.  APPLICATIONS 

Due  to  the  inefficiency  in  the  generation  of  nonlinear 
parametric  signals,  the  widespread  use  of  parametric  fog 
horns  is  not  practical.   It  does  provide  a  short  range  signal 
with  reduced  noise  pollution  that  can  convey  more  precise 
directional  information  and  with  scanning  can  warn  the  mariner 
of  danger  sectors  to  avoid  or  other  geographical  information. 
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